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Abstract 

In this short paper we show that making turbulence two- rather than tiiree- 
dimensional may increase effective turbulent viscosity by about 40%. Dimension- 
less hydrodynamical viscosity parameter up to otmax ~ 0.25A^i may be obtained 
in this approach, that is in better agreement with the observational data on non- 
stationary accretion than the values obtained in numerical simulations. However, 
the a-parameter values known from observations are still several times higher. 



1 Introduction 

Though a prescription jShakurj[T973l :[ Shakura & Svunvaev|[T973h proved to be 
quite useful for accretion disc physics, the physical mechanisms driving angular 
momentum transport are still not completely understood. Hydrodynamic as well as 
MHD turbulence and magnetic fields are generally proposed as the main viscosity 



sources. Une oi the sources or MML) turbulence is m agne to-rotationat ms t ability 
first predicted in the classic works bv lVelikhovl (Il959h and 'Chandr asekha3 h96(]|) 



and applied for the case of disc accretion by BeJbus & Haw lev ( 199^ 



Stationary accretion disc properties depend weakly on the viscosity a- 
parameter itself. Te// oc R~^^'^ effective t emperature profile ( and hence the out- 
coming multicolor blackbody spectrum, see iLvnden-BellH 19691) ) is independent of 
viscosity. However, most of the physical quantities (such as density, radial velocity 
and the total mass of the disc) depend somehow on the viscosity parameter. 

Detailed spectral fitting is capable for viscosity parameter estimates through 
the density of the disc atmosphere. Few attempts (such as . Davis et al..c2006,) ) were 
made yet to determine a-parameter from the X-ray spectra, but these estimates are 
still uncertain and highly model-dependent. 

Due to that reason, Q-parameter is generally determined from observational 
data on non-s tationary accretion phenomena such as dwarf nova fl ares and X- 
ray transients jSmakll 19991 : iKing et al.ll2007l : ISuleimanov et al.ll2008h . The actual 
viscosity parameter value generally affects characteristic viscosity timescales and 
subsequently the shapes of optical and X-ray lightcurves. lOiao & Lid J2009l) use 
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transitions between the low-hard and high-soft states in X-ray binaries to estimate 
the actual a-parameter values in the range 0.1 -j- 0.9. 

The dimensionless viscosity parameter is proven to vary in the range 0.1 -i- 
0.9 for hot, fully ionized accretion discs. Numerical MHD simulations generally 
predict a < 0.02 (see King et al. (2007) for review) much smaller than the values 
required by observations. 

Slightly higher val ues such as .1 -r- 0.2 appear in hydrodynamical turbulence 
models such as that o f lKatd jl994h . Turbulence is treated as a three-dimensional 
phenomenon. There are however reasons to consider also the two-dimensional 
case bec ause vertical turbulen t motions are damped by positive vertical entropy 
gradient i Shakura et alj 19781 : Tavler 19 8^). Here we follow the general line o f 
iKatd ( ll994l) : lNaravan et all ( [r994,) ; ,Kato & Yoshizawj ( ll997h ; lKato et al.l fl99^ . 
but consider two-dimensional case instead of three-dimensional. In the following 
section we write out the general form of the Reynolds stress evolution equations 
for the 2D-case. In section [3] a steady-state axisymmetric solution is established. 
We estimate the dimensional viscosity parameter value in section|4]and discuss our 
results and make conclusions in sections|5]and[6l 



2 Reynolds Stress 

Here we define the turbulent stress tensor as tij — {muj), 
fluctuating velocity components. The flow is considered incompressible, that 
presumably does not affe ct a non-divergent flow. We modify the equations by 
iKato & Yoshizawal ( Il997h in a way to account for the different dimensionality 
of the case considered. In our case, anisotropy tensor is naturally defined as 
bij = (tij — K)/2K, where z, j — r,ip and K = {trr + tipip)/2. Dissipation rate 
e should also appear without the 2/3 miltiplier present in the three-dimensional 
approach. Stress tensor evolution equations may be written as follows: 

d d \ 

+ ^t;— trr = 'iQtr^ + Urr — E (1) 



dt dip 
d d \ 

TTT ^ 1 — 2f2i(j5tp irr ^ Hy-tn (3) 

at dip J 2il 
For the pressure-strain correlatio n terms Hj, we adopt the general form with 
two dimensionless constants, used bv lK^ ( ll994h .n.. = -CiK>chj+C2KS^j. 
Ci,2 are universal dimensionless constants, and Sg = (dUi/dx j + dUj /dxi) /2 
is the ra te of strain of the avaraged flow. Following lKatd ( Il994l) . we use x rather 
than Q, ( Kato et al. 1998) as the characteristic frequency for local processes such 
as inertial waves. That assumption makes no difference for the case of Keplerian 
rotation when x = SI. Results may be modified by substituting CiQ/x instead 
of Ci in all the equations. We ingor e the two additional terms with coefficients 
Cs and C4 jKato & Yoshizawalll997l) . Their inclusion (for any C4 > and < 
C3 < 2) may only decrease the trip stress component value. 
Finally, pressure-strain correlation terms take the form: 

Urr = -^>c{trr-K) (4) 
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n^v = — (*^^ - K) (5) 
_ Ci C2 dO. 

— —Ktrip H ^fi^ (6) 

^ 2 ar 

3 Steady-State Axisymmetric Solution 

Stationary solution is easily obtained from equations l lU3t by zeroing their left- 
hand sides. It is instructive also to obtain an expression for the steady-case dissi- 
pation term as a sum of equations l|TJ and Kinetic energy K is also conserved 
in stationary case, hence: 

= - = -.-W-. (7) 

The dissipation rate may be therefore written as e = —rdQ/drtr^. Note that 
the two pressure-strain correlation terms exactly zero after summation. After some 
algebra one obtains the components of the Reynolds stress tensor (normalised by 
2K) as follows: 



/r, I 1 d 111 n \ dlnH N 
{"^ ' 2 dlnr J dlnr J 



2K 2 y Ao \ 2 1 



X 

d 111 n 



Ci n ( ^ C2 \ dluQ. 



2K AAq Q.\ 2 1 dlnr 



(8) 

(9) 

(10) 



where: 



4 Estimating l\irbulent a-Parameter 

Let us estimate the turbulent viscosity parameter for the two-dimensional case. We 
define the dimensionless a parameter as: 

p\tr^\ 7|irip| 

a = = — 5—, 

where Cs is the speed of sound, and Mt ~ 2K/cl is therefore the turbulent Mach 
number squared. For the general case of arbitrary Ci,2: 



7C1 X 





dlnO, 




dlnr 



Ml (12) 



As in the three-dimensional case, a-parameter has a maximum at a certain 
value of {Ci)max- In the two-dimensional case {Ci)max = 5 for Keplerian rota- 
tion. The maximal value of a is: 
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7 1 

^max J — 

O " /ill f!^ /dlllf2^2 

Maximal a-parameter values for two- and three-dimensional case differ by 
about 40%. In the particular case of Keplerian rotation, C2 = and 7 — 5/3, 

_ / 0.18A4? for 3D 
- I 0.257W? for 2D 

The second expression is the precise upper limit for two-dimensional hydro- 
dynamic turbulence. The difference with the three-dimensional case is that all the 
turbulent energy is concentrated in radial and azimuthal motions that contribute 
to the radial angular momentum transfer. Presence of the additional z dimension 
leads to decrease of radial and azimuthal motion amplitudes. 



dlnfi 



dliir 



Mi 



(13) 



5 Discussion 

The viscous a-parameter may be defined in different ways. Above we used a — 
Up = ptrip/p- Alternatively, a parameter is sometimes defined using the turbulent 
viscosity coefficient u = aiyC^/Q,. These two definitions of the a-parameter are 
related by the following expression: 



ptr^ _ 'yv 
P ci 







dlnO, 




dr 


= 7 


dlnr 


a„ 



(15) 



The two parameter definitions Up and a^ are absolutely equivalent. For Ke- 
plerian rotation and 7 = 5/3, a^ is higher by a factor 2.5. Therfore, our 
{cep),„ax = 0.25 corresponds to {a^),„ax = 0.1. 

Vertical structure of a thin accretion disc is unfavourable for strong turbulence 
development in z-direction. If net entropy increases with the distance from the 
equatorial plane of the disc, vertical motions are expected to be damp ed. That is the 
case for standard accretion discs in gas pressure-dominated regime ( Shak ura et al.l 
I1978| ; |t^y1 er 1980), excluding the zones of p artial ionization of hydrogen wher e 
super-adiabatic temperature gradient appears iMever & Mever-Hofmeiste3ll982h . 
Stable vertical stratification affects mostly the largest turbulent spatial and ve- 
locity scales, where vertical motions are effectively damped by buoyancy forces. 
The situation is different for horizontal turbulent motions that always have an en- 
ergy influx from the averaged shear flow. Thin accretion discs are predicted to 
be unstable with respect to radial convective motions (R o ssby wave instability, 
Morozov & Hoperskovl h990l) : iLovelace et all ( [199^) : Izhuravlev & Shakur j 



Largest vortices of the turbulent cascade receive their energy and mo- 
mentum immediately from the averaged flow. Due to that they are expected to be 
close to complanar with the accretion disc rotation. On the other hand, vertical 
motions do not interact directly with the mean flow. Because the mean flow is 
primarily affected by the largest eddies of the turbulent cascade, turbulent velocity 
field is effectively two- rather than three-dimensional, when its dynamical effect 
on the mean flow is considered. 

Excitation mechanisms and energy balance of hydrodynamical turbulent mo- 
tions in accretion discs are tightly connected with the instabilities present in the 
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flow. Though thin accretion discs generally fulfil the Reynolds and H0iland sta- 
bility criteria, they may be unsta ble to non-axisymmetr ic perturbations in the pres- 
ence of radial entropy gradient jLovelace et alj[l999h . Besides this, even in the 
linearly stable case there are certain ways to destabilize differentially rotating 
flows thro ugh non-linear instabilitie s, mode coupling and transient growth solu- 
tions (see lUmurhan & RegevI ( |2004) for review and references therein). Recur- 
rent transient growth solutions may lead to large turbulent energy growth for th e 
large Reynolds numbers characteristic for accretion discs jAfshordi et alJllOO^ . 
In short, it seems that turbulence is likely to develop even in linearly stable cases. 
Simulati ons are not consisten t however in this conclusion - certain numerical sim- 
ulations jHawlev et al.lll999h report stabilization of all the instability modes in 
shearing box simulations. 

Turbulence excitation in differentially rotating fluids are being extensively 
studied in laboratory experiments with Taylor-Cou ette flows. A comprehensive 
explanation of th e first experiments of that kind jTavloJ [l936h was given by 
IZeldovicri h98lh . who also points out the similarity between rotating fluid mo- 
tions and astrophysical discs. So far, laboratory experiments do not support high 
viscosity values in differe ntially rotatin g fluids with quasi-Keplerian profiles. Ef- 
fective a values found bv iJi et allfeOOdl) are well below ~ 10^^ for high Reynolds 
numbers Re > 10^. Taylor-Couette experiment, however, differs from astrophys- 
ical discs in some important points like radial entropy gradient and gravitational 
field of the central body. Resulting non-axisymmetric convective instability modes 
may be therefore responsible for high accretion rates in astrophysical discs. Be- 
sides this, astrophysical discs are highly supersonic flows, while Mach numbers in 
laboratory experiments are much lower unity. 

Numerical simulations mostly confirm the existence of long-living turbulence 
in rotating shear flows. Large scale two-dimensional vortices, mainly of anticy- 
clonic nature, ar e reported to form and survive (Bodo et al. 2007; Lithwick 2009|; 
IShen et al]|2006h . s ometimes exceedin g disc thickness in their spatial extent in az- 
imuthal direction. iBodo et sl\ j2007l) report that the azimuthal size of a typical 
vortex exceeds the scaleheight of the disc for low speed of sound values (or high 
Mach numbers of the flow), starting from about 0.05 in Keplerian velocity units. 
So the presence of the large-scale two-dimensional vortices relevant for angular 
momentum transfer is well confirmed by numerical simulations. It is not clear 
whether the vortices themselv es are stab le if the effect of the third (vertical) di- 
rection is taken into account. iLithwiclj 12009) claims that vortex dynamics and 
stability properties do not change significantly in 3D wit h respect to t he 2D case. 
On the other hand, another recent paper bv Lesur & Papaloizoul ( |2009|) report that 
elliptic 2D vortices are almost always destroyed by resonance processes with the 
local epicyclic oscillations. 

6 Conclusions 

Our main conclusion is that the turbulent viscosity Q-parameter may increase by up 
to 40% if the turbulence is two- rather than three-dimensional. The resulting upper 
limit {ap)max = 0.25Mt is still too low to explain the high viscosity parameter 
values derived from non-stationary accretion phenomena. 

It is possible that the contemporary turbulence theory based either on labora- 
tory experiments or a quasi-particle approach misses some important effects lead- 
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ing to more effective turbulent motion interaction. Alternatively, one may allow 
turbulence to be supersonic or to consider the influence of magnetic fields. The 
impact of the latter in controversial: as it was already mentioned in Introduction, 
MHD simulations predict relatively low a-parameter values. 

We conclude that the observational viscosity parameter is still 2-^3 times 
higher than the upper Umits given by theory. The problem is yet to be solved. 

This work was supported by the RFBR grant 09-02-00032. We thank the ref- 
eree (Marek Abramowicz) for useful comments and suggestions. 
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